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Abstract 

The electrostatic interaction between an excited atom and a diatomic ground state molecule in an 
arbitrary rovibrational level at large mutual separations is investigated with a general second-order 
perturbation theory, in the perspective of modeling the photoassociation between cold atoms and 
molecules. We find that the combination of quadrupole-quadrupole and van der Waals interactions 
compete with the rotational energy of the dimer, limiting the range of validity of the perturbative 
approach to distances larger than 100 a.u.. Numerical results are given for the long-range interac- 
tion between Cs and Cs2, showing that the photoassociation is probably efficient whatever the Cs2 
rotational energy. 
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I. INTRODUCTION 

In a recent paper (hereafter referred to as paper I) [l| we investigated the electrostatic 
interaction between an atom and a diatomic molecule at large mutual separations. This kind 
of study is particularly relevant in the context of the amazing development of researches on 
ultracold quantum gases, i.e. when the kinetic energy of the relative motion of the particles 
inside the gas is equivalent to a temperature far below 1 Kelvin. Many enlightening review 
papers are available in the scientific literature, like the most recent ones [2|-|4| devoted to 
ultracold molecular gases which is our main purpose. Motivations and applications of such 
researches are thoroughly discussed in these reviews, and extensive bibliography can be found 
there that we will not repeat in the present paper. Briefly, dilute atomic and molecular gases 
at ultracold temperatures exhibit pure quantum effects, as their dynamics is sensitive to 
quantum resonances and quantum statistics leading to Bose or Fermi degeneracy. Ultracold 
gases exhibit unique physical conditions for high precision measurements related to the limit 
of validity of fundamental theories like the Standard Model or for quantum simulation of 
condensed matter phases like superfluidity or superconductivity. 

At such a low energy regime, the relative motion of the particles inside the gas is controlled 
by their weak mutual interactions at large distance R, primarily induced by electrostatic 
forces described by a multipolar expansion in R~ n terms. For instance, it is well known 
that spin-free or rotation-less particles in their ground state interact mainly through van der 
Waals potentials varying as Cq/R 6 resulting from a second-order perturbation treatment 
of the dipole-dipole interaction. Calculations of the long-range dispersion coefficients Cq 



re species of choice for ultracold 
5|, |6[ and represent an invaluable 



for alkali-metal or alkaline-earth atoms -which represent t 
gases- have nowadays reached an unprecedented accuracy 
input for interpreting experiments. First-order terms may become dominant when one atom 
or both atoms are excited, resulting into interactions varying as R~ 3 Q or R~ 5 for 
dipole-dipole and quadrupole-quadrupole interactions, respectively. 

In principle, such calculations can be easily extended to the interaction between atoms 
and molecules with internal rotation. However most previously published studies have been 



restricted to the situation where molecules are fixed in space |9|, llOj. A recent study actually 



treated the van der Waals interaction between two identical ground state molecules in a 



given rotational level 



ll| using the second-order perturbation theory. In paper I, we were 



interested in the interaction between an excited atom and a ground state molecule in a 
rotational level, governed by a first-order quadrupole-quadrupole term varying as R~ 5 , that 
we evaluated through first-order degenerate perturbation theory. We demonstrated that 
this term competes with the rotational energy of the dimer, so that avoided crossings, and 
possibly long-range wells, could be expected in the long-range potential curves of the atom- 
dimer complex. Such patterns are relevant in the perspective of future studies aiming at 
associating an ultracold atom-molecule pair using laser light (photoassociation, or PA) to 
create stable ultracold triatomic molecules, according for instance to the reaction 

Cs(6P) + Cs 2 (X 1 S+) -»■ Cs^ -»■ Cs 3 + photon . (1) 

In the present work, we extend our treatment to the second order of the perturbation 
development, to determine the next term varying as R~ & which is expected to compete 
with the R~ 5 term when R decreases. As in paper I, we illustrate our development for the 
case of an alkali-metal ground state molecule in an arbitrary rotational level j with 

an alkali-metal atom excited to the n 2 P electronic level, but it can be easily generalized 
to arbitrary species. In Section [TT] we first recall the general expression for the long-range 
multipolar expansion and the expression of the Cq coefficients in the case of the van der 
Waals interaction between an atom and a molecule, and we relate these quantities to the 
dipole polarizabilities in imaginary frequencies, which are evaluated in Section UTTl Potential 
curves for the long-range interaction between an excited Cesium atom and a ground state Cs2 
molecule are displayed in Section IIVI before concluding remarks and perspectives (Section 
IV|) . Atomic units (a.u.) for distances (1 a.u. = 0.0529177 nm) and for energies (1 a.u. = 
219474.63137 cm" 1 ) will be used throughout the paper, except otherwise stated. 



II. EXPRESSION OF THE VAN DER WAALS INTERACTION 



We first recall briefly the notation used in paper I. The electrostatic potential energy 
between two charge distributions A (the dimer) and B (the atom) at large distance R (i.e. 



beyond the LeRoy radius 



121 ] ) is expressed as the usual multipolar expansion 
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where Q^ x {rx) is the operator associated with the 2 Lx -pole of charge distribution X (X = A 
or B), expressed in the coordinate system with the origin at the center of mass of X 



Qtirx) = ^ aZTTl £ «*f*y£ (k A) , (3) 

where 
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and L< is the minimum of La and Lb- The Z quantization axis for the projections ±M 
in the above equations is oriented from A to B, yielding the (— l) Ls factor in Eq. (j3J). We 
define two body-fixed (BF) coordinate systems (CS) displayed in Fig. [TJ the dimer CS (or 
D-CS) with axes X A , Ya, Za, and the trimer CS (or T-CS) with axes X, Y, Z. The T-CS 
is deduced from the D-CS by a rotation with an angle 6 around the Y axis. The T-CS is 
related to the space-fixed (SF) coordinate system (xyz) by the usual Euler angles (a,{3,j). 

The dimer is assumed in an arbitrary vibrational \vd) and rotational state \j) of its 
electronic ground state |X 1 S^"). The atom B is chosen with a single outer electron being 
excited to the p state \n,£ = 1). The projections rrij and A of j and I are defined with 
respect to the Z axis. The energy origin corresponds to an infinite separation between the 
atom and the dimer. The first-order perturbation theory yields the quadrupole-quadrupole 
energies with zeroth-order wave functions that can be written as a linear combination of the 
various \rrij,X) substates of the complex 



H) = ^2c mjX \m j ,X) (5) 

m,jX 

where rrij + A is a constant, and c mj \ are real. The Cq or van der Waals coefficient comes 
from dipole-dipole interaction whose expression V^{B) is obtained by setting La = Lb = 1 
in Eq. (EJ) 

VabW- e? ^ (1 + M)!(l-M)f {b) 

As the permanent dipole of both the atom and the dimer are zero for any states of the atom 
and the molecule, the (non-degenerate) second-order perturbation theory is used to obtain 





FIG. 1. The dimer [X^YaZa] (D-CS), and [XYZ] (T-CS) coordinate systems, defined for the 
dimer and for the trimer, respectively. The Za axis is along the dimer axis, while Z is oriented 
from the center of mass of the dimer towards the atom B. The Y and Ya axes coincide and point 
into the plane of the figure. The T-CS is deduced from the D-CS by a rotation with an angle 5 
around the Y axis. The Euler angles (a, /?, 7) are not represented here. 
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For a given |$q), the corresponding Cg coefficient can be expanded as a linear combination 
of crossed terms 

C 6 = C m Jl Ai £ V»y a A a C r 6 r > ( 8 ) 
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where the crossed coefficients Cg r are 
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The summation is performed over all possible states $° and $° of the A and B systems. 



With the approach of imaginary frequencies dipole polarizabilities (see e.g. Ref. 13|). 
the sum in Eq. (Q can be factorized into separated contributions from A and B. For this 
purpose, we use the identity 
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to transform the first term of Eq. (Q, with \x\ = E® — E° Xvd - and \y 
approach is applicable if E® — E\ Vd ^ > and E® — E® e > 0. This is always the case for the 
dimer if it is in its ground rovibronic state. Moreover, if the dimer is homonuclear, radiative 
transitions between rovibrational levels belonging to the same electronic state are forbidden, 
so that the identity still holds for any rovibrational level of the state. In contrast, 

atom B is in its first electronically excited state. Therefore, taking cesium as an example, 
Eq. ([TO]) is correct for all transitions, except Q 2 P — > 6 2 S, for which E^ — E® e < 0. In the 
latter case, a similar factorization can be performed, by setting \y\ = —E® + E® e , and by 
using the following identity 



du- 



2 \x\ 



\x\ — \y\ \x\ + \y\ x 2 — y 2 



xy 



+ 



7T J [x 2 + u 2 ) (y 2 + cu 2 ) 

2x 



x 2 — y 2 



(11) 



Inserting Eqs. f lTOj) and f lTTj) into Eq. (j9]), and following Ref. [14J, we obtain 
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where AE® 



El 



E® e , and G(x) is the Heaviside function, equal to 1 for each down- 
wards transition and otherwise. The expression of dipole polarizabilities is generalized to 

arbitrary frequencies (real or imaginary) according to 

(E° 



a MM' \ Z ) 



E°x VdJ ) 



x (Xv d jm h \Q[ 



M 



E° V 



Ql 



M 1 



Xv d jm 



.12 



a 



Ai A2 
-M-M 



{El 



(E° b - E^Y - z 



-M\ $ 
b 



-v2 
M' 



$0 



Q? 



niX 



(13) 



(14) 



x {ni\i \Q^[ 

where z can be real or imaginary In Eq.( fT2l the first term inside the brackets is the well- 
known product of the polarizabilities at imaginary frequencies of the molecular \Xvdj) and 
atomic \n£) states (the corresponding labels are omitted for clarity sake). In the second term 
appears the polarizability of the dimer at the real frequency of each downwards transition 
of the atom. We note that if the dimer polarizability is negative and significantly large, the 
second term in the square brackets can become larger in magnitude than the first term and, 
therefore can turn the Cq sign to positive. 



III. CALCULATION OF POLARIZABILITIES 



1 . The polarizability of the dimer 



Molecular polarizabilities are most often calculated ignoring rotation, so that it is wise 
to separate it in the equations of Section [III We start by expressing the dipole matrix 



element ^XA 



v d jm 



Qi 



M 



of Eq. (Tl3l) between the states defined with respect to the 
T-CS, in terms of matrix elements between states defined in the D-CS. For clarity, we 
explicitly write the dimer quantum number A, which is the projection of the electronic 



angular momentum on the Za axis. Specifying as well all the relevant quantum numbers 
according to |$°) = \X'A'v' d j'm'j^, we have 



(XAvtjmj Of 



(i 
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(15) 



where the quantities q± are the components of the electric dipole moment defined with 
respect to the D-CS. The index /x is either for E — > E transitions (A = A' = 0), or ±1 for 
E — t- n transitions. In the T-CS the rotational wave function of the dimer is proportional to 
d J m . Q and to d J m±1 , for E and II states, respectively. Using the properties of the d^ rotation 



matrices 



15| . we find 



(j m j \ d M»\f m j) 
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(16) 



where the C 3 ^^ m , and C^,_ (l are Clebsch-Gordan coefficients. If in Eq. (fT3|) we neglect 
the rotational part of the th -order energy (only low rotational levels are relevant for the 
cold temperature domain), i.e. 



E° Xvd + B Vd j (j + 1) « E 



Xv d i 



(17) 



we identify the exact expressions of the parallel polarizability a\\ and of the perpendicular 
polarizability a± with respect to the Za axis, for /x = and ±1 respectively. Therefore, Eq. 
( fl3|) becomes 

e|^[K°.)V) + 2 ( c s») 2 ^w 
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where we used the identity [15] 
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26+ 1 C7a ~ Q ' 

in order to put all the primes in the upper indices. The dependence on the rotational state 
of the dimer is restricted to the coefficients of the vibronic polarizabilities of the dimer. 
We note that Eq. (fT8|) is valid for real, imaginary, and zero frequencies (i.e. for static 
polarizabilities). 
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2. The valence contribution to the atomic polarizability 



We perform a similar development for the atom, in order to separate the radial and 
angular parts of the matrix elements of the dipole moment. By writing explicitly |$°) = 
\n'l'\') and by using Eq. ( I19p . we rewrite the dipole matrix elements as 



(ni\ 
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where r n ^ n ^i = (n£\r\n'£') is the matrix element of the valence electron radial coordinate. 
Eq. (TH1) becomes 
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The index £' above takes values corresponding to dipole-allowed transitions, i.e. £' — 0, 2 in 
the present case of a Cs(6 2 P) atom, and n' for all the relevant atomic levels. 

The similarity between Eqs. f|T8|) and ()2ip confirms the equivalence in the formalism be- 
tween the rotational angular momentum of the dimer and the electronic orbital momentum of 
the atom, which makes the generalization to more complicated cases like molecule-molecule 
long-range interactions quite straightforward. However, it is not possible to express Eq. (j2ip 
as a function of the sole polarizability of the atomic n£ level. For instance, if i = 1, M' = M 
and A 2 = Ai, the dipole polarizability <Xmm( z ) °f the sublevel nPXi contains angular fac- 
tors which are different for P — > S and P — >■ D transitions In the usual definition of 
the polarizability, an average is made over all the sublevels Ai leading to the disappearance 
of the angular factors which is not the case here. Therefore, we introduce state-to-state 

polarizabilities a n e n >e> associated to the different nl — > n'i' transitions 

2 {El, v - E°J 
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so that Eq. (121]) becomes 
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The state-to-state polarizabilities obey the property Y^n'i' a nirii' ~ QW> with a n £ being the 
actual (isotropic) atomic polarizability of the level nl. This identity is only approximate as 
the effect of core electrons have been neglected so far. 



3. The core contribution to the atomic polarizability 

Following Ref. js] we assume that the contribution of the core electrons can be treated as 
an additional correction a c to the total polarizability a n e, independent of the atomic state 
\n£) 

OLnl = «n<W + «c • (24) 

n't' 

In our numerical calculations, a c is obtained in the following way. First, the contribution 
of the valence electron is evaluated using tabulated values of dipole moments for transitions 
from the 6 2 S state of cesium. The result is then compared to the polarizability of Ref. j^] 
which accounts for the ionic core contribution, and the difference between our result and 
that of Ref. j^] is then attributed to a c . 

It is not straightforward to see that the core polarizability brings a simple additive con- 
tribution to to the Cq coefficients. As inner shells have different angular quantum numbers, 
we first consider only the electrons of a given closed inner shell, e.g. Ad for cesium. A sum 
over all the closed shells will be taken at the end. We note n c , £ c , tTj the principal, orbital, 
azimuthal and spin quantum numbers of the i th electron of the shell. The indexes n c and £ c 
are identical for the 2 (2£ c + 1) electrons, whereas hi varies from — £ c to 
Since the dipole operator is mono-electronic, the ionic core brings to the C§ coefficient an 



c and oi = ±\. 



additional term Cg, which is the sum of all individual contributions 
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where the general expression for Cg clc (i) is 
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In Eq. ( 126|) . the atomic states are characterized by independent electrons 

1' 
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and similarly |$°-) = \n c £ c , —£ c , — |) . . . \n' c £' c , k' i7 <Ji) . . . \n' c £' c , £' c , |) \n£\), and the molecular 
states, |$%) and |$°), have the same quantum numbers as previously. The states \&° B0 ) 
and |$^} are antisymmetric with respect to the permutation of two electrons. However, for 
simplicity, they are written here as simple products of the individual electronic states, which 
does not modify the value of Cq c1c . In Eq. (l26|) . the energy E® is independent on the index 
i, since i only makes a distinction between different projection of orbital and spin angular 
momenta. 

Unlike the valence electron, as the core is not affected by the C5 coefficient, there is no 
crossed term in the calculation of C^ clc (i). As a consequence, the Wigner-Eckart theorem 
imposes m'^ + M = rrij and m'- + M' = rrij, as well as k[ — M = ki and k[ + M' = ki, hence the 
condition M' = M. So, the double summation of Eq. ( )26l) reduces to a single one over M. 
Similarly to Eq. ( fl2l . we can factorize Eq. ( 1261) with polarizability-like quantities. Because 
there is no allowed transitions from the last inner shell of the atom to lower-energy states, 
that factorization does not contain any additional term due to the excited state. It yields 
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where a™^ 3 * s given by Eq. ( ITS]) and a n ^ 1 _ M is similar to f[2"Tj) : 
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In Eq. (129]) . the sum over k and cr runs over all the core electrons %. The sum over the two 
different spin projections for each orbital gives factor of 2. The sum of the Clebsch-Gordan 
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coefficients yields — | — . Therefore, the angular part of the dipole moment disappears and 
oi c _ M _ M finally reads 
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As one can see from the above expression, it is not necessary to separate different series of 
transitions l c —> £' c , and we obtain a meaningful polarizability. Finally, we can extend this 
results to all inner shells that gives the core polarizability 
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and the core contribution Cg to the Cq coefficient 
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With the approximations above, we obtain an additional term which depends on the physical 
polarizabilities of the atomic core and the molecule, and in which the rotation of the dimer 
appears explicitly. As all the factors in Eq. (132]) are positive, Cg is negative and, thus, it 
makes the interaction between the atom and the dimer more attractive (or less repulsive). 



4- General expression for Cq 

We summarize all the results of the previous paragraphs. We recall that the Cg coefficients 
are calculated for each eigenvector of the quadrupole-quadrupole operator. Each eigenvector 
|$g) is a linear combination of states \rrij, A) (see Eq. dSJ)) with rrij + A being constant. For 
a given |$q), the general expression for C 6 is 
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(34) 



is the atomic transition dipole moment, Q(x) is Heaviside's function, and AE®, e , nl = E®, e , 



E° 



IV. LONG-RANGE POTENTIAL CURVES FOR THE CS-CS 2 COMPLEX 



The atomic polarizability of Cs(6 2 P) comes from the accurate calculations of transition 
dipole moments in Ref. lfj, averaged over the 6P1/2 and &P3/2 levels. The atomic transition 



energies are extracted from Ref. 
the summation is restricted to n' 



17| . and averaged over fine structure levels. In Eq. (j22l) . 
7 to 10 for e = and ri = 5 to 10 for £' = 2. We 



use a mixture of experimental and theoretical molecular data to compute the molecular 
polarizabilities at imaginary frequencies, which involve summation over all vibrational leve 



of all electronic states of CS2. The CS2 ground state potential curve is taken from Ref. 



and the one of the B l Ii u {Qs + Qp) state from Ref. [191] . The A 1 H+(6s + 6p) and b 3 U u (Gs + 6p) 



and their spin-orbit coupling are those fitted to reproduce the data of Ref. |20| . All the 
other electronic states come from quantum chemistry calculations performed in our group 
according to the method described in Ref. 2l| as well as transition dipole moments which 
will be presented in a separate publication. The dissociation energies of the atom-dimer 
complex are given by CS2 rotational energies, B j (j + 1) where the rotational constant for 
the v d = level of Cs 2 is B = 1.17314 x 10" 2 cm" 1 Q. 
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Energy (cm ) 

FIG. 2. Calculated dipole polarizabilities in imaginary frequencies (a) for Cs2(X 1 £+ ,vj, = 0) (with 
the parallel and perpendicular components, full lines), and (b) for Cs(6 2 P) (lower panel, full line). 
The contribution of the sole B 1 H u (6s + 6p) state to a±, and of the pair of 0^ A^Ti^iQs + 6p) 
and 6 3 II n (6s + 6p) states to an, are shown (dot-dashed lines), (a): The polarizability components 
including the core polarizability are drawn with dashed lines, (b) The magnitude of second term 
of EqJT2] (dashed line), and of the core polarizability (dot-dashed line) are also shown. 

The variation of the dipole polarizability at imaginary frequencies is displayed in Fig. |2j 
The B l U u (Qs + 6p) state lowest electronic state contributes to the valence part for about 99% 
to a±, as well as the pair of states A 1 T>+(6s + 6p) and b 3 U u (6s + 6p) coupled by spin-orbit 
interaction to yield a pair of so-called 0+ states. The core polarizability represents only a 
few percents of the total one, for the represented frequency domain. We numerically checked 
that the dissociation continua can be neglected. 

The Cq coefficients are first calculated for Cesium in its ground state Q 2 S interacting with 
Cs 2 X 1 E+(w (J = 0, j) molecule, for which the quadrupole-quadrupole interaction is zero. The 
aw quantities now represent the actual dipole polarizability. Moreover there is no crossed 
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TABLE I. The Cq coefficients of the Cs2pT 1 £+, = 0, j)+Cs(6 2 5) long-range interaction calcu- 
lated for j = to 4. In analogy to a diatomic molecule, the Cq are sorted by projections of the 
total orbital quantum number mj = rrij on the Z axis (note that A = 0), and the parity ± through 
the reflection symmetry with respect any plane containing this axis. 



term in the C§ coefficients and Eq. ( 1331) reduces to 
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M j'm\ 
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■ 1 ^IMjrrij 



dua{iuj) 

2 



a±{tu 



(35) 



2j' + 1 

C( ) )\\\{iu) + 2(C{ 1 ) Q ) 

where a{iu) is the atomic polarizability in the 6 2 S' state, including core contributions. It is 
worth noting that in the particular case j = 0, Eq. fl35l) can be written in a similar form as 
for two S atoms 



doja(ioj)a(ioj) 



(36) 



with a = (ay + 2a±)/3 the so-called isotropic polarizability of the dimer. 

The results of our calculations are given in Table [I] and in Figure [31 They are sorted by 
values of the projection mj of the total orbital momentum on the Z axis, and for mj = 0, 
the parity ± through the reflection symmetry with respect any plane containing this axis. 
In analogy to diatomic molecules, these symmetries are labeled E ± , LT, A... All the Cq 
coefficients are negative, reflecting an attraction between the atom and the dimer which is 
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Atom-dimer distance R (a.u.) 



FIG. 3. Long-range potential energy curves Boj(j + 1) + Cq/R® of S + symmetry, as functions of 
the atom-dimer distance R (in logarithmic scale), representing the interaction between Cs(6 2 5) 
and Cs2(X 1 S+,^rf = 0, j) for j = to 4. The related Cq coefficients are given in Table HI 

slightly less than twice the one for two interacting Cs(6 2 S') atoms (6840 a.u. Q, 23]). 
This can be understood from the simple form of Eq. (I36[) . considering that the static dipole 
polarizability of a ground state Cs atom and of a ground state Cs2 molecule are respectively 
^402 a.u. 22 j . ^707 a.u. j^j], 25]. In addition, no curve crossing are visible on Figure |3]as 
all Cq coefficients have close values. Therefore, the validity of the long-range expansion is 
limited at short distances by the overlap of the electronic clouds. With (rg S ) = 42 a.u. 



we estimate the value of the LeRoy radius to Rlr ~ 40 — 45 a.u.. 

The potential energy curves for the long-range interaction between Cs(6 2 P) and Cs2(X 1 £+, 
0, j) are displayed in Figs. |4j[5l and|6j For the £ symmetries (see Fig. [4]), the curves calcu- 
lated in Paper I (only containing the C5 term) are also shown. The related C5 (from paper 
I) and Cq coefficients are listed in Table [III 

Due to large and negative Cq values, most potential curves are attractive. When the 
C5 coefficients are positive, tiny potential barriers are visible. For instance, the highest 
barrier is found for IT symmetry with a height of about 0.1 cm -1 in the curve correlated to 
the j = 1 rotational state. In paper I, we found that the low-i? limit of our perturbation 
analysis R m = 102 a.u. was imposed by the crossings between curves with different values 
of j. We see on Fig. HJ that the typical position of these crossings is only slightly modified 
by the Cq contribution. However, it is worth noting that the crossing range coincides with 
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FIG. 4. Long-range potential energy curves -Boj(j + 1) + C^/ R 5 + Cq/ R e (full lines), and Boj(j + 
1) + C5/-R 5 (dashed lines) as functions of the atom-dimer distance R (in logarithmic scale), for the 
S + and S~ symmetries, representing the interaction between Cs(6 2 P) and Cs2(X 1 £+, = 0, j) 
for j = to 4. The related C5 and Cq coefficients are given in Table UU 



the one in which the second-order contribution competes with the first-order contribution. 
As a consequence, the non-adiabatic couplings emerging in the crossing region, evoked in 
paper I, should contain both first-order and second-order terms. The resulting couplings 
will mix states characterized by a given symmetry (S ± , LT, ...) and a given j, with states 
characterized by the same symmetry and by f — j ± 2. 

Note that for the S + symmetry and j = 1, one see a long-range potential well, due to 
a negative C 5 coefficient (-1674 a.u.) and a positive C 6 coefficient (51249 a.u.) (see FigH]). 
This well is 0.9-cm -1 deep, and its minimum is located at R = 37 a.u., slightly below the 
LeRoy radius (~ 45 a.u.). We can expect that this long-range well will indeed exist, but 
that it will be modified by non-adiabatic couplings discussed above and by the electronic 
exchange. Therefore, its precise characterization requires quantum-chemical calculations. 
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V. CONCLUSIONS 



Using a second-order perturbation approach, we have shown that the long-range interac- 
tion between an excited Cs(6 2 P) atom and a Cs2(X 1 E+,t>d = 0, j) molecule is significantly 
modified compared to the first-order quadrupole-quadrupole interaction. Most potential 
curves exhibit an attractive behavior, except in a few cases where either a low potential 
barrier or a long-range well is visible. The barriers could prevent the collision to occur in 
the ultracold domain. The validity of the present approach is limited to a range of distances 
well beyond the radius defined by the conventional LeRoy criterion, due to the low energy 
spacings between molecular rotational levels. This is actually the main difference compared 
to atom-atom long-range interaction, which completely changes the physical conditions for 
PA of an atom with a molecule. The next step in the theory is to include non-adiabatic 
couplings between curves related to different rotational levels, which is currently in progress 
in our group. They will induce avoided crossings which may generate a complex collisional 
dynamics. 
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FIG. 6. Same as Figure [5] for the for the <3? symmetry (upper panel), and for V (solid line) and H 
(dashed line) symmetries (lower panel). 



While the present theory is general, it is difficult to predict if this situation is similar 
for all alkali systems, either homonuclear or heteronuclear, as it strongly depends on the 
balance between the influence of the various parameters, namely the quadrupole moment 
of the atom and the dimer, and their dipole polarizability. As already stressed in paper 
I, other effects should be taken in account, like atomic spin-orbit interaction. The fine 
structure splitting in Cesium (554.1 cm" 1 ) is much larger than the rotational energy of Cs2, 
so that this effect it will not dramatically change our description. However, it is important 
to stress that the above treatment was developed in the basis of the LS coupling case. The 
appropriate inclusion of the hyperfine structure will add to the complexity of the potential 
curves of Figs. HI |5j and El even it will not modify the main conclusion of the present study: 
the photoassociation of a ground state alkali- metal dimer molecule with a ground 

state nS 1/2 alkali- metal atom is generally possible by exciting the dimer- atom system with 
a laser field red-detuned from the nSi/ 2 —> nPi/2,3/2 atomic transition. 

19 



ACKNOWLEDGMENTS 



We thank M. Aymar for kindly providing us with Cs2 potential curves and transition 
dipole moments. This work was done with the support of Triangle de la Physique under 
contract 2008-007T-QCCM (Quantum Control of Cold Molecules), and National Science 
Foundation under grant PHY-0855622. 



[1] M. Lepers, O. Dulieu, and V. Kokoouline, Phys. Rev. A 82, 042711 (2010). 

[2] L. D. Carr and J. Ye, New J. Phys. 11, 055009. 

[3] O. Dulieu and C. Gabbanini, Rep. Prog. Phys. 72, 086401 (2009). 

[4] C. Chin, R. Grimm, P. Julienne, and E. Tiesinga, Rev. Mod. Phys. 82, 1225 (2010). 

[5] J. Mitroy and M. W. J. Bromley, Phys. Rev. A 68, 052714 (2003). 

[6] A. Derevianko, S. G. Porsev, and J. F. Babb, At. Data Nucl. Data Tables 96, 323 (2010). 
[7] N. Bouloufa, A. Crubellier, and O. Dulieu, Phys. Scripta T113, 014014 (2009). 
[8] M. Marinescu and A. F. Starace, Phys. Rev. A 55, 2067 (1997). 

[9] B. Bussery-Honvault, F. Dayou, and A. Zanchet, J. Chem. Phys. 129, 234302 (2008). 
[10] B. Bussery-Honvault and F. Dayou, J. Phys. Chem. A 113, 14961 (2009). 
[11] S. Kotochigova, arXiv:1003. 2672(2010). 
[12] R. J. Le Roy, Can. J. Phys. 52, 246 (1974). 

[13] D. Spelsberg, T. Lorenz, and W. Meyer, J. Chem. Phys. 99, 7845 (1993). 

[14] C. Zhu, A. Dalgarno, S. G. Porsev, and A. Derevianko, Phys. Rev. A 70, 32722 (2004). 

[15] D. A. Varshalovich, A. N. Moskalev, and V. K. Khersonskii, Quantum theory of angular 
momentum (Leningrad, 1988). 

[16] U. S. E. Iskrenova-Tchoukova, M.S. Safronova, J. Comput. Meth. Sc. Engin. 7, 521 (2007). 

[17] C. Moore, Atomic energy levels, vol. 2 (Chromium through Niobium) (U.S. government print- 
ing office, Washington, 1952). 

[18] C. Amiot, O. Dulieu, R. F. Gutterres, and F. Masnou-Seeuws, Phys. Rev. A 66, 052506 (2002). 

[19] U. Diemer, R. Duchowicz, M. Ertel, E. Mehdizadeh, and W. Demtroder, Chem. Phys. Lett. 
164, 419 (1989). 

20 



[20] J. G. Danzl, E. Haller, M. Gustavsson, M. J. Mark, R. Hart, N. Bouloufa, O. Dulieu, H. Ritsch, 

and H.-C. Nagerl, Science 321, 1062 (2008). 
[21] M. Aymar and O. Dulieu, J. Chem. Phys. 122, 204302 (2005). 

[22] A. Derevianko, W. R. Johnson, M. S. Safronova, and J. F. Babb, Phys. Rev. Lett. 82, 3589 
(1999). 

[23] A. Derevianko, J. F. Babb, and A. Dalgarno, Phys. Rev. A 63, 052704 (2001). 
[24] J. Deiglmayr, M. Aymar, R. Wester, M. Weidemiiller, and O. Dulieu, J. Chem. Phys. 129, 
064309 (2008). 



[25] The quantity 2ce c (Cs + )=30.8 a.u. deduced from Ref. 
[26] M. Aymar, private communication. 



22f | has been added to the one of Ref. 



24]. 



21 



symmetry 


3 


C 5 (a.u.) 


C 6 (a.u.) 


symmetry 


3 


C 5 (a.u.) 


C 6 (a.u.) 


£+ 








-42704 


n 


4 


-739 


671 




1 


-1674 


51249 




4 


108 


-15884 




1 





-21562 




4 


522 


-47279 




2 


-913 


12128 


A 


1 


-279 


-18694 




2 


116 


-16885 




2 


-140 


-21244 




3 


-796 


4923 




2 


1136 


-95614 




3 


145 


-15420 




3 


-835 


-1624 




4 


-755 


2251 




3 


-87 


-19563 




4 


157 


-14835 




3 


736 


-65454 




1 





-43920 




4 


-721 


-2643 




2 


399 


-45131 




4 


-11 


-17153 




3 


465 


-45333 




4 


623 


-56200 




4 


489 


-45407 




2 


-399 


-16589 


n 








-23605 




3 


-245 


-18030 




1 





-29303 




3 


1175 


-103161 




1 


1116 


-79756 




4 


-783 


-5057 




2 


-964 


7305 




4 


-161 


-17444 




2 


-19 


-22961 




4 


835 


-76003 




2 


584 


-50736 


r 


3 


-465 


-15420 




3 


-783 


2496 




4 


-320 


-16392 




3 


64 


-17295 




4 


1208 


-107555 




3 


532 


-48103 


H 


4 


-507 


-14676 



TABLE II. The C 5 and C 6 coefficients of the Cs 2 (X 1 'E+,v d = 0, j)+Cs(6 2 P) long-range interaction 
calculated for j = to 4. In analogy to a diatomic molecule, the states are sorted by projections 
of the total orbital quantum number mj = rrij + A on the Z axis, and the ± reflection symmetry 
through any plane containing this axis. The values for C§ are taken from Paper I. 
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